Gradient marker" — a universal wave pattern in inhomogeneous continuum 
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Wave transport in a media with slow spatial gradient of its characteristics is found to exhibit a 
universal wave pattern ("gradient marker") in a vicinity of the maxima/minima of the gradient. 
The pattern is common for optics, quantum mechanics and any other propagation governed by the 
same wave equation. Derived analytically, it has an elegantly simple yet nontrivial profile found in 
perfect agreement with numerical simulations for specific examples. We also found resonant states 
in continuum in the case of quantum wells, and formulated criterium for their existence. 
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Wave patterns in inhomogeneous media or confining 
structures are of great interest to quantum mechanics, 
optics and electrodynamics, acoustics, hydrodynamics, 
and chemistry. Examples include wave packets in atoms 
[1], Ghladny patterns in acoustics, EM resonator and 
waveguide modes [2] , Anderson localization in disordered 
systems [3], soliton formation [4] due to nonlinearity, in- 
cluding atomic solitons in bosonic gas [5], as well as gi- 
ant waves near caustics [6], waves in chemical reactions 
[7], dark-soliton grids [8], "scars" in "quantum billiard" 
[9], "quantum carpets" in QM potentials [10], nano- 
stratification of local field in finite lattices [11], etc. In 
all of those, the presence of multi-modes or a broad-bend 
spectrum is pre-requisite for interference and pattern for- 
mation in inhomogeneous or confining structures. 

In this Letter we show, however, that a localized wave 
pattern - an immobile single-cycle intensity profile - can 
emerge in a single- mode wave in a vicinity of a min/max 
of the gradient of QM potential or optical refractive in- 
dex. The phenomenon is universal for both optics and 
quantum mechanics, and for any other propagation de- 
scribed by a wave equation (1) below. What makes it un- 
usual is that it emerges in media with no potential wells 
and only a smooth inhomogeneity yielding no reflection, 
- and is originated by a purely traveling wave with ap- 
parently no other modes to interfere with. We found, 
however, that this wave here generates a co-traveling but 
localized " satellite" of slightly different phase and ampli- 
tude resulting in " self-interference" . The wave ideally is 
not trapped and carries its momentum and energy flux 
unchanged through the area. To a degree, the pattern 
mimics a 2-nd order spatial derivative of the refractive 
index (or potential function); it would be natural to call 
it a "gradient (G) marker". In QM it may be most pro- 
nounced for an above-barrier propagation of electron in 
continuum over smoothly- varying potential; in solid state 
it might emerge above the critical temperature for the 
Anderson localization to vanish. Even for a potential 
well, when the energy of electrons exceeds the ionization 
potential and there is no trapping, the G-markers emerge 
as the main non-resonant localized feature. 

To demonstrate the effect and elucidate analytical re- 



sults (to be compared with numerical simulations) we 
consider ID-case written, for the sake of compactness, 
in "optical" terms, using space- varying refractive index 
n(x); yet we consistently "translate" all the effects and 
approaches into QM-terms. A ID spatial dynamics of 
an w-monochromatic plane wave with linearly polarized 
electrical field E = e p E{x) exp(—iuit) + c.c, propagating 
in the a;-axis (here e p _L e x is a polarization unity vector), 
is governed by wave equation 



E" + n 2 (0E = 0; £ = xk ; 



(1) 



where fco = oj/c = 27r/Arj, and "prime" stands for d/dt;. 
(For H field, e x _L H _L e p , one has H = —iE' in non- 
magnetic materials; for a traveling wave, \H\ = n\E\, 
if n = const.) In QM-terms, this corresponds to ID- 
scattering of a particle in continuum by a potential U (x), 
with E replaced by a wave function, ip, of a particle, H 
- by {-iyj') = p Q /k h, n - by p c /k h, where p Q and 
Pc = \/2m[Eo — U (x)} are its quantum and classical 
momenta respectively, and Eo - full energy. We will 
consider only the case n 2 > 0, where one can attain a 
no-reflection mode of the main interest to us here; oth- 
erwise, with n 2 (£) crossing zero, the system may exhibit 
a full reflection characterized by an Airy function as e. 
g. near a turning point in QM [12], or a critical point in 
plasma [2], or caustics in optics and water waves [6]. 

Eq (1) is ubiquitous in physics and engineering. Since 
few known functions n(£) allow for analytical solutions, 
numerical simulations and/or approximate analytical so- 
lutions in general have to be used. Of the most interest to 
us here will be the limit of adiabatically slow variation 
in space, when gradient parameter /i ~ (fcrji'imm) -1 , 
where L is a spatial scale of inhomogeneity, is small, 
fx <C 1, which corresponds to a quasi-classical case in QM. 
The reflectivity R in this case vanishes as R = 0{e~ A ^) 
[12,13], where A = 0(1) (usually A > 1), and reflection 
can be neglected by a large margin. A solution is pro- 
vided then by a WKB approximation [12] as traveling 
waves, C^' exp(±i J ndQ/n 1 / 2 . Considering e. g. a 
forward wave, and setting n' — > at |£| — > oo, where we 
normalize its intensity by setting fil-El 2 ^ = 1, we look 
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for next approximation as a perturbed WKB solution 

E= [1 + A(0]e i J" n */n 1 / 2 with A = 7 + i/3 (2) 

where A(£) (|A| 2 <C 1) is a slow- varying complex func- 
tion, A — > at |£| — > oo, 7 and ft - real; as we will see 
later on, \j\ max ~ |/3'|max = 0(^i 2 ). Using ansatz (2) 
in Eq. (1), setting real and imaginary parts of the sum 
of all the perturbations terms to zero, and collecting the 
terms of lower order in \x in each one of them, we obtain 
for the real part an equation consisting of 0(/x 2 ) terms 



$ = -(n'/n 3 / 2 )'/^ 1 



/2 



(3) 



and for the imaginary part - an equation consisting of 
0(/x 3 ) terms, integration of which yields 7 = — (f3 2 + 
(3'/n)/2, where we set the integration constant to zero 
due to above condition A± OQ = 0. (It is worth noting 
that in the end, all the terms with 1 get canceled, so 
there is no need for further integration of Eq. (3).) We 
can finally arrive at a G-marker intensity by calculating 
the perturbation, 51 (£) = 7 — 1, of the normalized field 
intensity I = n\E\ 2 = (1 + 7) 2 + /3 2 w 1 + 2 7 + /3 2 , 
retaining the terms lowest in /x, and obtaining to o(xt 2 ): 



Sift) = -p'/n = (n'/n^ 2 )'/4n 



3/2 



(4) 



In the vicinity of a gradient peak, <57(£) makes an asym- 
metric single-cycle shape, with its middle point shifted by 
0(koL) toward the area with lower refractive index (or 
higher potential); its higher (and positive at that) peak 
is also located in the same area, see Figs. 1 and 2. One 
can see that 8I(£) more or less mimics a second deriva- 
tive of n. Eq.(4) can also be obtained via quasi-classical 
approximation in QM [12], whereby one has to search 
for high-order corrections for the phase of ip as function 
of the classical momentum pc, after which it has to be 
translated into correction to intensity. 

How far the asymptotic result (4) can be pushed be- 
yond the limit /iCl, and what is a critical /x cr = 0(1), 
can be explored only by numerical simulations, which 
also helps to reveal a real nature of a small parame- 
ter zx (which appears to be substantially different from 



a standard (|n'|/n 



<C 1 [14]). Before comparing 



Eq. (4) to numerical simulations for specific profiles n(£) 
and various /x, let us make sure it conforms to the conser- 
vation of EM energy flux, i. e. (time-averaged) magni- 
tude, S, of the Poynting vector, S — E x 77/2 (in Abra- 
ham's form) in general case. Writing S = (Ee~ lut +c.c.)- 
(77e _l "* + c.c.)/2, and ^-averaging it, which amounts here 
to omitting terms with e ±21 "* , we have S = Re(EH*) 
= Re(iE*'E). In QM terms, it corresponds to mathemat- 
ical expectation of a particle momentum, < VIpqIV' >; 
Pq = —ihd/dx. Using Eq. (2) and retaining the terms of 
the lowest (2-nd here) order in zx, we have S = I + (3'/n; 
due to Eq. (3) it confirms that S = 1 = inv to o(/x 2 ). 

For numerical simulations of Eq. (1) with an arbitrary 
profile n(£) and arbitrary /x, we broke it into two 1st order 
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FIG. 1: (a) Refractive index, n (and potential U) soft-step 
spatial profiles, m = 1.5 and ri2 = 3, fj, C r ~ 0.24; (b-d) G- 
marker intensity, SI, vs distance x/2L for various parameters 
H; curves: 5In - numerical, and SI a - analytical, Eq. (4). 
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(Maxwell, in EM-case) equations: E' — iH; H' = in 2 E. 
To model a "soft step" n(£), we use a function (Fig. la): 

n(x) = m + (n 2 - m)[l + tanh(2a;/i)]/2 (5) 

with controllable L, n,, and its gradient parameter as 

^ = .g(n 1 - 1 +n 2 - 1 ); g = (AqL)" 1 = A /2^L (6) 

[14]; the spectral dispersion of rij can be safely ignored 
here. The calculations with an arbitrary [i are done by 
using a multi-point algorithm and a "reverse propaga- 
tion" mode, whereby we start at £ 3> ,g _1 , postulate 
that only one (transmitted) field remains there, E^ — > 
exp(jri2^)/v / ^2' (^e Sommerfeld's condition), and then 
go backward, till reaching a symmetrically located area 
in front of the gradient, £ <C — where we record the 
intensity of an incident wave, !,„ = ni\E+H/rii\ 2 /4, and 
then normalize all the stored intensities by 7j„ [15]. The 
precision of calculations is checked against the deviation 
of S, at each point from that recorded at the incidence; 
typically it was better than 10~ 6 . 

The retroreflection from the gradient area is strongly 
suppressed and G-markcr is well emphasized (see e. g. 
Figs, lc and Id) provided that /i < ji cr , where parameter 
fi cr was found by us to be almost universal, [i cr s=s (27r) _1 
at r n = ni/n 2 + n 2 /ni ^> 1, and slightly increasing to 
ji cr 1/4 near n\ ~ n 2 . The highest (positive) G- marker 
peaks, SI max , can easily reach a few percent of the inten- 
sity, I, especially at r n 3> 1. In the case of a "shallow" 
soft step, \n\ — 7T-2 1 <C n/2, n = (n\ + n 2 )/2 [in QM this 
would correspond to a kinetic energy Eo much higher 
than the drop of potential, Uo, E /C/o ~ n/\n\ — n 2 | 3> 1], 
the max/min of 81 are almost of the same magnitude, 

\SI M \^2g 2 \ ni ~n 2 \/(5 3 / 2 n 4 ) (7) 

and located at x w ±L/4. (In general, the parallel be- 
tween optics and QM can be guided by the relationship 
L/ /E = l-min[(nf/ni),(n|/n?)].) 

Figs. lb-Id for the case of m = 1.5, n 2 = 3, [i cr — 
0.24, show numerical simulations of spatial dynamics of 
51 n, converging amazingly fast to an asymptotic ana- 
lytical result for a G-marker intensity, SI a, Eq. (4), as 
soon as \i < [i cr . Fig. lb with fi/fi cr = 3/2 shows a 
residual reflection giving rise to an oscillating structure 
(partial standing wave), comparable in its amplitude to 
a G-marker, while Fig. lb (ji — fj, cr ) depicts distinct and 
strong G-marker formed even at L w A /2. Finally, an 
inhomogeneity with L = A (Fig. Id, fi/fi cr = 2/3) is 
sufficient to produce a very clean G-marker. 

We move now to investigate G-marker formation by a 
potential well (or a refractive index plateau) by modeling 
it with an "up-and-down" double step, Fig. 2a: 

n(x) = m + (n 2 - n 1 )(T+ - T_)/2 tanh(D/L) (8) 

where T± = tanh[(2a; ± D)/L], and D is a controllable 
spacing between the steps. For D <C L, it becomes 



n(x) = rti + (n 2 — n\)/ cosh 2 (2a;/L), but for our purposes 
here we choose more box-like well, D/L = 8, which has 
/i defined by Eq. (6), and the same pL cr as a soft-step (5). 
As expected, both walls form G-markers symmetric to 
each other, Figs, lb and lc, so that to form a G-marker 
it does not matter which way a wave is arriving - from the 
lower index or from the higher one. At /i > fj, cr one can 
see some oscillations, same as for a single wall in Fig. lb 
for the same /i, and ideally clean G-markcrs for /j, < fi cr , 
similar to Fig. Id for the same /i. 

The major difference here comes, however, in the area 
fj, > fi C r- Here, at certain (countable) set of points in the 
continuum, while there are strong oscillations within a 
potential well, which indicates a significant wave reflec- 
tion between G-markers, there is no reflection from the 
entire potential well, see Fig. 2d. That confined par- 
tially standing wave is a signature of a resonant state in 
a finite-depth quantum well with rigid walls, most known 
in the case of a finite rectangular box. Fig. 2d depicts 
one of those states with E /?7o = 4/3. The condition 
for them to emerge above a quantum well is a significant 
rigidity of the well's walls, /x > fx cr . In the limit /i 3> /i C r, 
their energies in the continuum coincide with those of a 
finite box, or in turn - with the eigenstates of a box with 
infinitely- high walls, En = (Nhn) 2 /2mD 2 , where N is 
a natural number, provided that En > Uo. In optics 
terms, they correspond to full-transmission resonances of 
a Fabri-Pierrot resonator with semi-transparent mirrors. 
In solid-state, these states may reveal themselves dur- 
ing a (5-kick field ionization via production of spatially- 
stratified bunches in photoelcctron current, whose kinetic 
energies coincide with those of the resonant states [16]. 

Potential uses/applications of ID (or almost ID) G- 
markers can be envisioned, such as (a) observation of 
quantum "traces" in continuum, i. e. beyond "quantum 
carpets" [9] in potential wells, (b) detection and con- 
trol of slight changes of optical fiber parameters [17], (c) 
the diagnostics of cold under-dense plasma, (d) medical 
surface- wave ultrasound tomography, (e) detection of the 
movement of near-shelf profiles of the bottom of oceans 
and rivers by space- or air-borne photography of the pat- 
terns of wind-driven gravitation waves, as well as (f) 
contour-detection and tracing of submerged large/long 
moving man-made objects or whales in the ocean. 

A 2D and 3D expansion of the theory may need to be 
developed for other potential applications of G-markcrs 
such (g) the "tomography" of quantum landscape in dis- 
ordered solid-state at above-critical temperature [3], (i) 
a bulk tomography of opaque fluids (e. g. oil or muddy 
water) by using non-penetrating surface EM or acoustic 
waves, or of solid-state bodies (e. g. in " introvision" of 
computer chips, or lacunas in blobs of metallic alloys or 
glass), as well as (j) in plasma- and astro-physics. 

In conclusion, we predicted the formation of a univer- 
sal feature in wave transport in a inhomogeneous media 
- a standing single-cycle spatial modulation of wave in- 
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FIG. 2: (a) Refractive index plateau, n (and potential well 
U) profiles, with m, 712, and /^ cr same as in Fig. 1; (b,c) (57, 
vs x/2L for /i > /i cr (b), and /i < fi cr (c), with and SI a 
as in Fig.l; (d) resonant state in continuum (see the text). 



tensity - gradient marker - located in closed vicinity of 
max/min of a gradient of refractive index or potential 
function. We found a critical condition for a G-marker 
to be resolved on the background of residual reflection. 
In the presence of a trapping potential, we also found 
resonant states/modes in the continuum at the energies 
above photo-ionization and formulated the condition of 
those modes to exist when G-markers are not dominant. 

This work is supported by AFOSR. 
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